Fluctuations of the Casimir-like force between two membrane inclusions 
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Although Casimir forces are inseparable from their fluctuations, little is known about these fluc- 
tuations in soft matter systems. We use the membrane stress tensor to study the fluctuations of 
the membrane-mediated Casimir-like force. This method enables us to recover the Casimir force 
between two inclusions and to calculate its variance. We show that the Casimir force is dominated 
by its fluctuations. Furthermore, when the distance d between the inclusions is decreased from 
infinity, the variance of the Casimir force decreases as —1/d?. This distance dependence shares a 
common physical origin with the Casimir force itself. 
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In 1948, Casimir predicted that two uncharged metallic 
plates placed in vacuum should attract each other [1]. In- 
deed, the boundary conditions imposed by the plates on 
the electromagnetic field constrain its quantum fluctua- 
tion modes in such a way that the zero-point energy of the 
system depends on the distance between the plates. Such 
long-range fluctuation-induced forces arise between any 
objects that impose boundary conditions on a fluctuating 
fleld with long-range correlations [2] . Casimir-like effects 
driven by thermal fluctuations of material media were 
first discussed by Fisher and de Gennes in the context of 
critical mixtures [3-5]. They also appear in soft matter 
systems such as liquid crystals [6, 7], fluid membranes 
[8-12] and fluid interfaces [13], as well as in superfluids 
[14]. The first direct measurement of thermal Casimir- 
like forces was achieved very recently in a critical binary 
mixture by Hertlein et al. [15]. Experimental evidence 
for fiuctuation-induced forces between lipid domains in 
vesicles was also provided very recently [16]. 

The Casimir force actually coincides with the average 
of the stress exerted by the fiuctuating field [17]. Al- 
though Casimir forces are by essence inseparable from 
their fluctuations, the latter have been scarcely studied. 
In 1991, Barton first characterized the fiuctuations of the 
quantum Casimir force by calculating the variance of the 
stress tensor [18] . There have been few studies since then 
[19]. The fluctuations of Casimir forces are, however, 
of fundamental importance. Indeed, Casimir force mea- 
surements are always performed by probing a fluctuating 
quantity, cither the force itself [20, 21] or the position of 
one of the interacting objects [15]. In addition, the dis- 
tance dependence of the fluctuations of Casimir forces is 
intriguing as it shares a common origin with the Casimir 
effect: the suppression of fluctuating degrees of freedom. 

The study of the fluctuations of Casimir-like forces was 
initiated by Bartolo et al., who considered the case of par- 
allel plates imposing Dirichlct boundary conditions on a 
fluctuating scalar field [22]. Although the case studied 
in Ref. [22] is quite generic, including, e.g., classical spin 
systems [23], it does not encompass all the soft matter 
systems that can give rise to a Casimir-like effect. Mem- 



brane inclusions, for instance, which are small objects 
interacting with the membrane through a complicated 
stress tensor [24, 25], arc left out. 

In this paper, we calculate the fluctuations of the 
membrane-mediated Casimir-like force acting between 
two inclusions, e.g., proteins, that locally constrain the 
curvature of the membrane. First, we show that the 
Casimir-like force, which is usually derived from the 
membrane free energy [8~11], can be obtained by inte- 
grating the average membrane stress tensor. Then we 
calculate the variance of this force and we study its dis- 
tance dependence. We also discuss the effect on the vari- 
ance of the interplay between the Casimir-like force and 
the curvature-dependent force. 

Model. — Biological membranes are fluid lipid bilayers 
that are well described on a large scale by an elastic cur- 
vature energy known as the Helfrich Hamiltonian [26]: 



H= /dr|[V^Mr)]^ 



(1) 



where h{v), with r = {x,y), denotes the height of the 
membrane with respect to a reference plane, and k is the 
membrane bending rigidity. This form holds for symmet- 
ric membranes, constituted of two identical monolayers, 
undergoing small deformations with respect to the flat 
shape and for distances smaller than ^k/ct, where a is 
the membrane tension. The components of the stress 
tensor S associated with the effective Hamiltonian H 
are given by [25]: 



E - - 

- 2 



{dpiY^idlhY +K{d^h)d^v^h, (2) 



(3) 



We consider inclusions that locally constrain the cur- 
vature of the membrane. This corresponds to the 
most generic case for inclusions not subject to exter- 
nal forces or torques. Particles such as proteins with 
Bin/amphiphysin/Rvs (BAR) domains [27] and viral 
capsids [28] are well-known examples of inclusions pro- 
ducing local membrane curvature [29-31]. 
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FIG. 1. (Color online). Atomistic molecular dynamics simu- 
lation of a BAR domain inducing a strong local curvature in a 
lipid membrane. Reproduced with permission from Ref. [27]. 
Copyright (2006) National Academy of Sciences, USA. 



FIG. 2. (Color online). Two curvature-inducing inclusions of 
radius a separated by a distance d in a fluctuating membrane. 
The black line is the contour used for the calculation of the 
Casimir force by integration of the stress tensor. 



We model the inclusions as point-like objects. This is 
justified since the typical radius of membrane proteins 
is comparable to the membrane thickness (see Fig. 1), 
which vanishes in our coarse-grained description. This 
simplification makes the calculation of the force variance 
tractable. A disadvantage of this model, however, is that 
the size of the inclusions and the ultraviolet cutoff are 
not independent from one another. 

Average force. — We first consider the case of two in- 
clusions imposing a vanishing curvature, since it gives 
the Casimir-like force without the curvature-dependent 
contribution that would arise from an average deforma- 
tion of the membrane [8, 10, 11]. In order to calculate 
the average of the stress tensor, we need the correla- 
tion function of the membrane height, and therefore the 
partition function. The latter is obtained by integrat- 
ing over all the configurations satisfying the constraints 
d^h = dxdyh = dyh = at the positions ri and y-2, of the 
two inclusions: 



Z[u] = Vh Y[S {B^h{R,)) e 



-l3H+fdr h{r)u{r) 



(4) 



i=l 



Here, /3 ~ and u{y) is an external field conjugate 

to /i(r); we have used the vectors D ~ (9^, dxdy, dy, d^, 
dxdy, dy) and R = (ri, ri, ri, r2, r2, r-i) to express the 
constraints. Writing the Dirac distributions as Fourier 
transforms and integrating first over h{v) and then over 
the Fourier variables yields [10, 11]: 



Z[u] = Z[0] exp 



i j drdr'u(r)C(r, r')M(r') 



(5) 



where C(r,r') = {h{v)h{r')) = G(r - r') + G"(r, r') 
is the correlation function of the membrane height. 
G(r) — /q s*'* is the Green function of the op- 
erator associated with the Helfrich Hamiltonian (1). 
Here and in the following, we use the shorthand = 
{k-QT/K) Jj^^^^i^i^^ (iq/(27r)^, where A is the ultraviolet 
cutoff, comparable with the inverse membrane thickness, 
and ^ stands for the membrane size in the absence of 



tension, or the fluctuation correlation length ^ « a/k/ct 
otherwise. The second term in the correlation function, 
which comes from the presence of the inclusions, is found 
to be 

6 

G'(r,r') = - ^ [D,G(r - R,)] M^T^ D,G(r' -R,), (6) 



where Mij ^ D,DjG(R, - Rj). 

Let / be the force exerted by the rest of the membrane 
on a circular membrane patch of radius r centered on the 
inclusion number 2 (see Fig. 2). Its component along the 
X axis joining the two inclusions is 



(r) cos 9 + Y;xy{Y) sin 



(7) 



with r = r (cosfl, sinfl). Let us calculate the average 
of fx- Thanks to Eq. (2-3), {T^xx) and (Say) can be 
expressed as linear combinations of the derivatives of the 
correlation function G. For instance, the first term of Yixx 
gives {{dyh{v)Y) = dydy,C{r,r')\yi=y. The divergence of 
the average stress tensor vanishes everywhere except on 
the inclusions, which means that {fx) is independent of 
the contour chosen as long as it surrounds one inclusion. 
Indeed, for all r smaller than the distance d between the 
inclusions, we obtain the Casimir-like force 



(/,) = 24A:bT 



(2/A)^ 



o{{My 



(8) 



which was first derived from the membrane free energy 
[9-11]. This result coincides with the Casimir-like force 
between two rigid disks of radius a [8] for A = 2/a [9-11]. 
Thus, a = 2/A can be interpreted as the effective radius 
of our point inclusions. 

Fluctuations. — Using the stress tensor provides a way 
to calculate the variance of fx'. 

A/^ ~ da dd Vxx (r, r') cos 6 cos a 

Jv L 

-|-2]4y(r,r')cos0sina -|- V^y(r,r')sin0sinQ; , (9) 
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= r (cos a, sin a 



and V = 



where Fy(r,r') = (E„(r)E 
with r = r (cos 6*, sin 0), r' 
[0,27r] X [0,27r]. 

Given the partition function (5), we may use Wick's 
theorem to express the correlation functions Vij in terms 
of C. For instance, the term in 142;(r,r') that originates 
from the product of the first and second terms in E^.^, 
gives Vxx{Y,r') — [9y9|,C(r, r')] . We represent 

diagrammatically its contribution to A /J as 



r~ / da d9 Vxx (r, r') cos 6 cos a = 



(10) 



Each hne stands for the propagator C(r,r'), and the 
dashes (resp. dots) come from dx or dx' (resp. dy or 
dyi) depending on whether they stand on the left or on 
the right of the diagram; the angular integrations are 
understood. All the diagrams that contribute to A/J 
are displayed in the supplementary material [32]. Since 
C = G+G", each diagram breaks into four terms. For the 
diagram in Eq. (10), one of these terms reads explicitly: 



9^fc2A,(q)Af^lA,(k) q'^q'^ 



q.k.q' 



/ dad6'cos6lcosae*('i+i') ''e'(*'-i') ''', (11) 



where the double line stands for G'(r,r') while the 
dashed line stands for G(r — r'), and Aj(q) ~ QjC^^"^'^' 
with Q = {ql, qxqy, q^, ql, qxQy, q^)- The integrals 
in Eq. (11) can be performed analytically in the un- 
physical limit A?' <C 1. The result will guide our 
understanding of the physical case r = a. Summing 
all the diagrams in A/^ [32] then yields Af^ ~ 
r4A6(/fcBT)2 [In (A^) - 4{Ad)-^ + O {{Ad)-^)] /192,' 
where zero-average oscillations at the cutoff frequency 
with amplitude in {Ad)~^/^ have been discarded. 

Since even a piece of inclusion-free membrane is sub- 
ject to a fluctuating force of finite variance, A/J depends 
on the contour chosen to calculate fx, contrary to (fx)- 
Thus, in order to obtain the fluctuation of the Casimir- 
like force acting on an inclusion, we should take a contour 
that includes the inclusion and only it. As in our model, 
the effective radius of the inclusions is a = 2/ A, the best 
we can do is to choose r = a and to set, from now on, 
A = 2/a. In order to calculate A/J for r = a, we have 
computed the integrals of the diagrams numerically. Our 
study (see Fig. 3) shows that the leading behavior of A/J 
in a/d is well described by the formula 



A/2 - 0.112 



In 



2^ 



0.239 



(12) 



which has the same main features as the analytical re- 
sult found in the limit Ar <C 1 . The only difference is the 
presence of the constant —0.239, which is negligible com- 
pared to the logarithmic term. As it can be absorbed in 
the deflnition of ^, we shall discard it from now on. In the 




FIG. 3. Left: Non-dimensional force variance A/| = 
/3^a^A/^ for d — i- cx) as a function of ^/a. The points are 
results from our numerical integrations, the line is the fit 
A/f = 0.112[hi(25/a) - 0.239]. Right: d dependence of A/f 
for S,/a = 10"^. Oscillations at the cutoff frequency have been 
smoothed out. The line is the fit A/|(d — oo) — A/| = 
0.112 aVrf^ + 7.18 a'^/d'^. Similar results were obtained for 
other values of ^/a (10^ and 10*), with the same coefficient 
in front of jd^ . 



physical regime, corresponding to ^ >> a and d > 2a, 
(12) is always positive. 

Our result (12) shows that A/J is dominated by a 
distance-independent term, which corresponds to the 
fluctuations of the zero-average force exerted on a sin- 
gle inclusion by the membrane bulk. To leading order in 
a/d, the signal-to- noise ratio for the Casimir force is thus 



24 



0.1121n 



2C 



(13) 



so that the Casimir force is very small compared to its 
fluctuations in the physical case d > 2a. Besides, the dis- 
tance dependence of A/^ originates from the suppression 
of fluctuation modes by the boundary conditions imposed 
by the two inclusions, so it shares a common physical ori- 
gin with the Casimir effect. To extract this d dependence, 
we can define the Casimir effect relative to the fluctua- 
tions as 



dd 



0.224 



d3 



(14) 



This result, obtained for point inclusions, is independent 
of the cutoff. The question whether this universality 
holds for extended domains with dimension a indepen- 
dent of the cutoff A would be interesting to address. 

Inclusions inducing a non-zero curvature. — When the 
inclusions impose a finite local curvature (see, e.g.. 
Fig. 1), a component due to the average deformation 



4 



of the membrane adds to the Casimir force [8, 10, 
31]. The partition function can be obtained using the 
same method [10, 11] as in the zero-curvature case dis- 
cussed before. The correlation function {h{r) h{r')) — 
(/i(r))(/i(r')) remains unchanged, equal to C(r,r'). The 
only change is the average deformation of the membrane: 
(/i(r)) = Vi I)jG{r -Kj), where V contains the el- 
ements of the curvature tensors of the two inclusions. 
We restrict ourselves to inclusions imposing the same 
isotropic curvature c ^ a~^, so that V = (c, 0, c, c, 0, c). 
The correspondence with finite-size inclusions imposing 
a contact angle variation a over a length a is obtained 
for c = a/a. 

Integrating the average stress tensor over a contour 
that surrounds one inclusion but not the other (sec Fig. 2) 
gives the sum of the Casimir force (8) and a curvature- 
dependent force, i.e.. 



(/,) =24fcBT 



d5 



(15) 

We thus recover from the stress tensor the total 
membrane- mediated interaction [8-10]. 

Studying the fluctuations of is a straightforward 
generalization of the zero-curvature case. The only dif- 
ference is that Wick's theorem applies to h{r) — {h{r)) 
instead of h{r), so that diagrams involving (/i(r)) ap- 



pear in the calculation. In the limit Ar ^ 1, we 
obtain analytically A/J - r'^A'^{kBT)^ [(t>ia)ln{Ai) - 
4V'(a)(Ad)-2-HO((Ad)-4)]/i92 with 0(a) ^ l + l2nl3Ka^ 
and '4'{(y) = 1 -l- Stt/Sko^. The conclusions drawn in 
the vanishing curvature case still hold, but the Casimir 
effect relative to the fluctuations becomes dAf^/dd ^ 
h{Ar)'^{kBT)'^ip{a)d~^ . Thus, in the non-small curva- 
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ture regime a > {SttISk) we obtain 



~dd~^ 



d^ 



(16) 



The k-BTKo^ factor, which replaces the (fceT)^ appearing 
in the pure Casimir case, reveals the interplay between 
the Casimir force and the curvature-dependent force. 

In summary, the stress tensor is a powerful tool that 
allows for studying the Casimir-like force between mem- 
brane inclusions, especially its fluctuations. Using a 
coarse-grained description in which the inclusions are 
point-like, we have calculated the variance of the Casimir 
force. Our results show that the fluctuations dominate 
the average force, and that they depend on the distance 
between the inclusions. Further possible developments 
include treating the case of extended inclusions, and test- 
ing our results using coarse-grained membrane numerical 
simulations such as those in Ref. [31]. 
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